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We report the application of the nonlinear er model to study the multi-skyrmion problem in the 
quantum Hall ferromagnet system. We make use of a first-principle calculation to derive an ana- 
lytical form for the inter-skyrmionic interaction to show that the ground state of the system can be 
described by a ferromagnet triangular Skyrme lattice near v = 1 where skyrmions are extremely 
dilute and a continuous transition into antiferromagnet square lattice occurs by increasing the 
skyrmion density and therefore \v — 1|. Using these results we demonstrate that the transition for 
a triangular to a square lattice which was previously derived, using the Hartree-Fock method, can 
also be seen in the field theory picture. We investigate the possibility that the skyrmions bound in 
pair to make a bi-skyrmion triangular lattice when the Zeeman energy is extremely small. We show 
that the energy of a skyrmion with charge Q is less than the energy of Q skyrmions each with charge 
one when the short range interaction among them is considered. By taking the quantum fluctua- 
tions into account, we also argue the possibility of the existence of a superconductor-insulator and 
the non-zero temperature phase transitions. 
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I. INTRODUCTION 

Recently, theipe have been a number of experimentalEHil 
and theoreticaHia investigations of skyrmions in the 
integer and fractional quantum Hall effect (QHE) 
in 2-dimcjpsional electron gas (2DEG). It has been 
reportea£3c3 that the microscopic calculations show the 
ground state of a two-dimensional electron system at the 
Landau level filling factor near u-^ 1 is a Skyrme, (soli- 
ton) lattice. Recently Brey et al&B and Cote et alt3 have 
made use of the time-dependent random phase approxi- 
mation (RPA) to show that the zero-temperature ground 
state of a multi-skyrmion system is a centered square lat- 
tice. In particular they have found the possibility of the 
structural phase transition into triangular Wigner crys- 
tal by changing the Zeeman energy and the Landau level 
filling factor. One may expect on the general grounds, 
these microscopic observations should be supported by 
other many body techniques, e.g., the appropriate field 
theory. It has been shown by Sondhi et alE and Moon et 
al.U, that the quantum Hall ferromagnet can be described 
by a generalized minimal nonlinear a (NLer) model to 
add the Zeeman and Coulomb energy, as an alternative 
approach. The consistency of this minimal field theory 
with other microscopic calculations for a system of a sin- 
gle skscrmion has been reported-numerically by Abolfath 
et al.n Recently, Green et alt3 have exploited a varia- 
tional approach on the generalized NLcr-model to study 
the multi-skyrmion system. They have indicated that 
the ground state is closer to a hexagonal structure, i.e., 
a distorted centered square lattice. Surprisingly, they 
came to a quite different conclusion from those obtained 
by microscopic consideration. One naturally seeks the 
reasons for the reported discrepancies between the dif- 



ferent approaches employed for the multi-skyrmion sys- 
tems. The aim of this paper is to present a proper field 
theory picture which avoids such controvertial issues. By 
exploiting th|e.jpinimal field theory of the quantum Hall 
NLcr-modeH'BB we present a description for the field the- 
ory which is consistent with the microscopic picture. In 
this way, we present a first-principle calculation to derive 
the analytical form for the inter-skyrmionic interaction. 
To the best of our knowledge this type of calculation 
has not been reported before. In the past, several au- 
thors have used a phenomenological approach to guess 
the form of the interaction by using the known Hactree- 
Fock long wave length fRin waves spectrum .cJO'EII In 
contrast to Green et aZ.til, we explicitly show that the 
ground state is a centered square lattice and/or triangu- 
lar Wigner crystal, deponding on the value of the Zeeman 
energy and/or Landau level filling factor. These results 
which werp previously derived, using the Hartree-Fock 
method,E2rn3 can also be seen in our field theory picture. 
We start from a proper superposition rule for the unit 
vectors of the NLcr-model to evaluate the correct analyt- 
ical form for the inter-skyrmionic interaction. We there- 
fore demonstrate that the field theoretic description of 
the NLcr-model does indeed lead to conclusion in agree- 
ment with microscopic Hartree-Fock calculations. We 
take advantage of the particle-hole symmetry^ (which is 
equivalent to skyrmion-antiskyrmion duality in the clas- 
sical field theory) to study the skyrmions associated with 
holes (y < 1). We will show that within v\ c < v ~ 1, 
the ground state is a triangular lattice in order to re- 
duce the Coulomb energy, since the effective interaction 
energy between skyrmions through the channel of the 
Zeeman term is negligible. Increasing the skyrmion den- 
sity in such a way that << v^ c < v < v\ c , leads to the 
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formation of a centered square lattice. The square lattice 
is favored by the effective interaction made by the gradi- 
ent and Zeeman energy. The effective interaction which 
is the attractive XY model allows the antiferromagnetic 
alignment of the skyrmion's orientation. This can be 
anticipated by the frustration of the skyrmionic hedge- 
hog fields within the triangular lattice which is the favor- 
able orientation for the Coulomb energy. This leads to 
a structural phase transition due to varying the Landau 
level filling factor. The attractive XY interaction makes 
a tendency for collapsing the single skyrmions which are 
located in different positions. It favors recombination of 
the iV-seperated single skyrmions and forming a single 
skyrmion with topological charge N. The competition 
between the attractive XY interaction and the repulsive 
Coulomb interaction leads to stability of the Skyrme lat- 
tice of the single-skyrmions. We will discuss for any short 
range interaction, like the hard-core model, the strength 
of the repulsive interaction is not sufficient to make a 
stable lattice. In this case, the energy of a charge-Q 
skyrmion is lower than any other configuration, hence the 
destruction of the lattice. We demonstrate in special sit- 
uations, at extremely small values of the Zeeman energy, 
(depends on the filling factor) a single skyrmion square 
lattice (with antiferromagnets orientation) collapses to a 
triangular lattice of bi-skyrmions (charge-2 skyrmions). 
In the latter, the orientation of skyrmions is governed by 
the ferromagnets ordering. By taking the quantum fluc- 
tuations into account, we also argue the possibility of the 
existence of a superconductor-insulator and a non-zero 
temperature Kosterlitz-Thouless phase transition, Tkt- 
We also show that the 2DES with the typical Zeeman en- 
ergy at filling factor smaller than one (the validity range 
of the square lattice) are in the insulating phase. We fi- 
nally estimate the melting temperature, T m , for the lat- 
tice via the Kosterlitz-Thouless-Halperin-Nelson- Young 
(KTHNY) mechanism. A comparison between Tkt and 
T m shows that both mechanisms can be accounted for 
these systems, depending on the Zeeman splitting factor 
and/or the Landau level filling factor. 



II. AN OVERVIEW ON NLcr-MODEL IN QHF 

In the NLcr- model approach to QHF, the spin of an 
electron may be described classically by a unit vector 
(spin coherent state), as an order parameter, whose di- 
rection may be changed continuously in the space. In this 
representation, the effective Hamiltonian is functional of 
the unit vector m(r)E 



E[m] = E [m] + E z [m] + E c [m], 



(1) 



where Eq [m] and E z [m] are the conventional NLcr model 
and the Zeeman energy respectively. The last term i? c [m] 
is the long range Coulomb energy due to the equivale: 
between Hall current and skyrmion topological curren 
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2n£ 2 



d 2 r [ 1 - m z (r) 



E c [m] = - 



d 2 r 



iV p(r)p(r') 



(2a) 



(2b) 



(2c) 



The long range nature of the Coulomb interaction is cru- 
cial for stability of th e Sky rme lattice. We will get back 
to this point in Sec. |III B| . Here p s = e 2 /(16v / 27re^o) is 
the spin stiffness at v = 1 (assuming zero layer thick- 
ness for the 2DEG) which arises from the Coulomb ex- 
change energy, e is the background dielectric constant 
of the semiconductor, g = ge 2 /(2e£o) is the Zeeman 
term, g is the effective gyromagnetic ratio, and £orf& 
the magnetic length. The charge density is given byOEl 
p(r) = -g^r e Q /3 m(r) • [<9 Q m(r) x dptn(r)] which is equal 
to the filling factor times the topological 0(3) spin tex- 
ture density of the quantum Hall ferromagnet (QHF). 
The total skyrmion charge denoted by Q is an integer- 
valued topological invariant. It can be determined by the 
integration upon the charge density and classified by ho- 
motopy group of a 2D-sphere respectively. The lowest en- 
ergy skyrmion solution, m(r), has to satisfy a non- linear 
differential equation which can be obtained by minimiz- 
ing the energy in Eq.(nh with respect to m 



p s V 2 + m ■ V 2 m) m M - -^—^ (5 zft - m z m 
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e al3 {d a V(r)}(m x dpih)^ = 0, 



(3) 



where V(r) is the Hartree potential (the exchange poten- 
tial resides in p s ) 



e J r-r' 



(4) 



and p is the skyrmion charge density associated with the 
minimum energy solution, m(r). The solutions of Eq.(||) 
can be classified by the skyrmion charge Q = J cPr p(r). 
From now and for the sake of simplicity, we remove the 
tilde over the classical solution and denote it by m. It 
is easy to find the following equation of motion of the 
optimal skyrmion by making use of the cross product of 
m upon Eq.(||) 



d a J a = x m )A, 



(5) 



where 



J a = p s (m x d a m)\ - —V{Y)e a pdpm\. (6) 

47T 

One may immediately read off from Eq.(||) that d a J„ = 
0, i.e., a conserved iso-spin current. We may define the 
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ground state of QHF at v = 1 as a vacuum of skyrmionic 
spin textures where all spins are aligned along the mag- 
netic field direction, i.e., the z-axis. Note that in the 
absence of the Zeeman energy, the alignment of spins 
along an arbitrary axis occurs due to spontaneous global 
0(3) symmetry breakinga hence the minimizing the elec- 
trons exchange energy. This is also the state of spins far 
from the center of the skyrmions. Note that the num- 
ber of skyrmions (antiskyrmions) in the ground state is 
counted by \N — N$\. Before we consider a lattice of 
skyrmions we need to have the correct shape of a sin- 
gle skyrmion. A single skyrmion is a topological optimal 
solution of Eq.(|3[) with Q = 1. For our purposes, it is 
convenient to parameterize the unit vector m by 

m = {<p x ,<Py, \Jl-ipip), (7) 

where ip — ip x + itp y . Near the core of the skyrmion we do 
not expect that the shape of the skyrmion is influenced 
much by the magnetic field. However, for large distances 
from the core, the Zeeman energy becomes dominant and 
we need to consider its effect where the direction of the 
unit vector m is close to the vacuum, namely, z-axis. 
Taking the limit of small ip, we can expand the m up to 
quadratic order in the ip. Putting this in Eq.(f?]) gives 

2#] = J d 2 r (^VV + 4^2^) , (8) 
which leads to the equation of motion 

-, s V^ + ^ = 0. (9) 

Introducing k 2 — g / (2it£qP s ) we simply have the equa- 
tion — V 2 ip + K 2 ijj = 0. We are interested in the 'vor- 
tex' solution ip = 2d z T, where d z = (d x + id y )/2 and 
z = x + iy. (If g = 0, this would result in V 2 T = with 
the solution T oc ln(r) and therefore ip oc z/r 2 ). Substi- 
tuting this we find -V 2 T + k 2 T = or T oc e'^/y/r 
and therefore that ip oc ze~ Kr /r 3 / 2 far from the core. 
The most important of this part is of course the exx>p=. 
nential (as opposed to, algebraic if g = 0), fall-off BE£ 
To check this, we solve Eq.(||) for a single-skyrmionr 
The results obtained from these numerical calculation 
shows 0(r) = \ip(r)\ = cKKi(nr) /2ir far from the core 
of skyrmion where K\(x) = \J-k / (2x)e~ x (1 + 0(l/x)) 
is the modified bessel function, c is a constant, can be 
obtained from the asymptotic form of 9(r) 

c = 30.4. (10) 

The dynamics of a skyrmion spin texture in NLa-modcl 
may be incorporated via the Wess-Zumino action.ES 
The result of expansion for the single skyrmion's Wess- 
Zumino term is 

Swz = 2 J drj dim- (d T m x 9 t m) 

= 4f <*r/A«r, T )^(r, T ), (11) 



where we keep the quadratic terms in Eq.(pT|). At this 
level of approximation, the effective action may be ob- 
tained by Eq. (||) and Eq. ( J]~T| ) where the Lagrangian den- 
sity is 

and S = J dr J d 2 rC. The optimal solution of ip is 
then identical to the solution of the time dependent 
Schrodinger equation where the external potential is pro- 
portional to the Zeeman splitting factor. Then the single 
skyrmion behaves approximately like a quantum mechan- 
ical point particle far from its core. Eq. (|l2|) describes the 
usual gapful ferromagnetic spin wave mode. 

III. SKYRMION INTERACTION 

We can now consider the interaction between 
skyrmions by generalizing our linearized energy func- 
tional. For that one should note that for the optimal 
single skyrmion spin texture in Eq. (^) we may choose a 
particular orientation. In general the texture can be ro- 
tated without costing any energy. Since the system shows 
the U(l) symmetry, any valid skyrmion spin texture can 
be obtained by rotating all spins about the £-axis by an- 
gle x- Therefore the state, m' = exp(ix/ z )m, is also an 
optimal solution of the Hamiltonian where I z = —id/dx 
is the generator of the rotation along the £-axis in the in- 
ternal space. One may expect that it contributes to the 
Hamiltonian through (h 2 /2A )(— id/dx — Co) 2 where A 
is the moment of inertia of the single skyrmion. This 
is the leading term of the total energy which is ex- 
panded with respect to the number of reversed spin, £. 
In general the dimensional analysis shows E z oc £ and 
E c oc The optimal value of the number of reversed 

spin, Co = E z /(2g), can be evaluated by minimization of 
the total energy with respect to £ 

Zo=^2 J d 2 r[l-m z (v)}, (13) 

where m z (r) is the optimal solution of Eq.(||) corre- 
sponding to the given Zeeman splitting factor. This 
leads to the optimal Coulomb and Zeeman energy and 
then the predicition of E c /E z = 2 and hence Ao = 

h 2 (d 2 E/d£ 2 )^ = h 2 E z /{6g 2 ). The validity of this sim- 
ple dimensional analysis has been demonstrated numeri- 
caly in Ref. 9. Note that £ is a quantized variable since 
X is compact. 

A. Square lattice 

Our goal is now to calculate the interaction energy be- 
tween skyrmions with different orientations. Following 
Piette et a/.li3, we start with the conventional NLcr model 
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to find out the proper superposition rule for the hedge- 
hog fields of different skyrmions. In the absence of the 
Zeeman and Coulomb energies, the energy functional is 
scale invariant and one may find the optimal solutions an- 
alytically. In this case, it can be shown that any-analytic 
complex polynomial defines an optimal solutiontJilj 



m(r) 



2w x 



2w„ 



1 



w ■ 



1 



\w 



2 ' 1 



tv 



(14) 



where w = w x + iw y is any Q-sector analytical function. 
One may decompose w into a series of analytical func- 
tions each with Q = 1 



JV 



(15) 



Any Q = 1 analytical function represent a single- 
skyrmion, then the number of skyrmions N is clearly the 
total skyrmions topological charge, i.e., N = Q. Eq.(|l5|) 
denotes a sequence of order parameters nx, (r) in configu- 
ration space. It is convenient to parameterize the nij(r) 
by 

mj(r) = (smr}j(r) cos Q(r),smr]j(r) sia(j(r), cost] j(r)) , 

(16) 

where r\j (r) and Q (r) are polar and azimuthal field 
variable— associated with the jth skyrmion. One may 
defineaO Q = ip — Xj where ip is the standard azimuthal 
angle, e.g., r = (r cos ip, r sin ip) and Xj is skyrmions ori- 
entation and measuring the deviation from the standard 
hedge-hog fields. In Fig. [l] the in-plane relative orien- 
tation, Xi between two skyrmions with relative distance, 
R, is shown. If we have a gas of skyrmions far away from 
each other, the total energy is invariant under variation of 
skyrmions orientation. For finite separation we expect a 
coupling between skyrmions due to different orientations. 
Here we consider a situation where are localized and 
well separated. This is valid for dilute skyrmions in a 
quantum Hall system, e.g., v-i c < v < v\ c - One may 
divide the configuration space (R 2 ) into N regions such 
that Ui is significant in region i and small in others 



m 



1 + z ■ m 4 



(17) 



m = m i + !l 1 x m; + -Q* x x m») + 0(fi 3 ), (20) 

here Eq.(|20|) describes an infinitesimal rotation of 
about Qi axis where 

fk = im, x {(1 + z ■ mO^ff - (m, • & cS )z}, (21) 



and 



N 



** ff (r)=X>ito- 



(22) 



Therefore the effect of the other skyrmions on the specific 
skyrmion is the same as a single skyrmion with charge 
Q — 1 due to the effective linear field, <i> c ff. In order 
to study the effect of the Zeeman and Coulomb energy, 
we make the ansatz that the above superposition rule 
is valid for skyrmions far from each other even in the 
presence of the full interaction. This can be taken into 
account by Eq.(pO|) and evaluating the total energy of 
multi-skyrmion in QHF. It is obvious that the total en- 
ergy can be divided into energies in separated regions. 
We assume that the interaction between the skyrmions 
is weak, hence 17^ = ir^ x 3>g ff . One may obtain easily 
the total energy by redoing the same calculation for all 
separated regions and sum over energies 



E[m] = 



N 

E 

i=l 



E[mi] 



9 



Anil 



cir 



(23) 



where ip\ s = <& l x cS + i® l y off is a complex field and £^[1x1] 
is the effective Coulomb interaction between skyrmions. 
Note that in the absence of the Coulomb term, the sad- 
dle point solution associated with the scalar field, ip* , are 
vortices, i.e., — V 2 ?^ -I- K 2 tp l = 0. One may divide the 
total energy, Eq.(E3T), into two parts, the self energy of 
skyrmions and the interaction energy which are desig- 
nated by T and V respectively 



T[m] 



N 



£[m, 



d 2 v 



N 



9 



4< 



(24) 



Here /i = (a;, y) and 



where j ^ i and $j = (tp J x , <pP , 0) 



'j 1 > 



1 + v/1 - Wj ■ 



(18) 



(19) 



Substituting Eq. (p^|)-Eq. (pj|) into Eq.(14) and expanding 
nij up to tpj, we may find m in the ith region 



The first term in Eq.(gJ) is the total self energy of the 
isolated skyrmions and the second term is the effect of 
their tail in other regions, i.e., the contribution of their 
energy in regions far from the core. Here we are inter- 
ested to study the effect of interaction in a system of 
many skyrmions and their physical relevance. One may 
find the effective interaction between the skyrmions by 
making use of the Stokes theorem and neglecting the 
next nearest neighbor terms in Eq.(|23|) which are de- 
scribed by the terms like ^(Pr^j^^J^k^j ' ® k anc ^ 
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V cS {m]=E° oS {m]+Eg[ m ] 1 



(25) 



and 



£ e ° ff [m]=E Urd a {J^n iX } 
i=i Ji 



(26) 



<ii> 



where jW has been defined for the ith skyrmion and 
satisfies the Euler-Lagrange differential equations, Eq.(j|) 
and Eq.(^). ££ ff [m] is the contribution of the gradient 
and Zeeman energy to the effective interaction. It de- 
scribes a system of interacting dipoles (see Fig. ||). Ex- 
panding the charge density of skyrmions, p, in terms of 
ip, leads to the effective Coulomb interaction 



E^[m]^^Jd 2 rJd 2 r 



l Pi {v)p {v') 
Ir-r'l 
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+ [d 2 r d^V l {r)d y m J • <>, • ...} . (27) 

where Pi = j^f e a fl • [<9 Q m.; x dpmi\. The first and sec- 
ond term in Eq.(|27|) are the Coulomb energy due to the 
monopole and dipole counterparts of skyrmions respec- 
tively. The former falls off like Rr 1 whereas the latter 
falls as R~ 2 where the distance between two skyrmions is 
denoted by R. Unlike the dense skyrmions in which the 
effect of the dipole term is crucial, the monopole term 
dominates for an extremely dilute (i.e., v ~ 1) skyrmions 
system. Neglecting the Coulombic dipole term, one may 
evaluate the integrals in Eq.(p6|) by applying the tech- 
niques thatp-were developed for a pair of skyrmions by 
Piette et aUB 

2. [ j2„/Pi( r )Pi( r ') 



V e n[m]= e -Y J jd 2 rjd 



2 ~ 

+ H E cos fe - Xi)*o(K|Ri ~ R il)> (28) 



2 ~ 

<_9_ 

4tt 2 

<v> 

where Xj ~ Xi an d Rj — Ri describe the in-plane rela- 
tive orientation and effective distance between skyrmion 
the jth and the ith. The first term in Eq.(|2|) is the 
electrostatic monopole term independent of the relative 
orientation. The normalization factor, c = 30.4, is calcu- 
lated for a single skyrmion numerically. The second term 
in Eq. ( p8| ) describes a classical XY-model where the min- 
imum energy configuration specifies the relative orienta- 
tion corresponds to \j ~ Xi — 7r - Ko(%) is the modified 
Bessel function, i.e., the coupling between the site i and 
j decays exponentially. We see that there is an exponen- 
tial decrease of the coupling between the ith and the jth 



skyrmions for R 1/k. We will show that an antiferro- 
magnet ordering of the hedge-hog fields within centered 
square lattice is the minima of the effective potential, 
Eq. (^8|)E£l, in some circumstances although the Coulomb 
interaction in Eq. (pq) is long range and the XY term con- 
tributes to the total energy as a short range interaction. 
The stability of the square lattice can be anticipated via 
generalizing the standard techniques for calculating the 
collective modes pjihich has been developed by the Bonsall 
and MaradudineO for the multi-skyrmion system] 



B. Recombination of the skyrmions 

One may consider the static and dynamical properties 
of a skyrmion lattice using the interaction that we have 
derived in the previous section 



v; ff = E^>(|Ri-Rj|) 



+ 2 J(|Ri-R,-|)coB(xi-X7)- (29) 

<ij> 

In our model Vo(R) is the electrostatic monopole inter- 
action and J{R) = j^Kq(kR). In Fig. ||, J is depicted 
as a function of the Landau level filling factor for a given 
Zeeman energy. The sign of V e s specifies the global min- 
ima of the energy functional. For instance, if V c g is pos- 
itive then E(Q = N) > NE{Q = 1) and the system can 
be composed by seperated single skyrmions. The valid- 
ity of this inequality has been investigated recently by 
Lilliehook et alEB when N > 2. The case of N — 2 is 
marginal and will be discussed in the next subsection. 
The ionization energy of the lattice per skyrmion is iden- 
tified by E(Q = 1) where R — > oo. Here E(Q) is the 
optimized energy of an individual Q-skyrmion (the self 
energy of the skyrmion). One may immediately notice 
that the XY term in Eq. (^9|) favors forming a Q-skyrmion 
(Q = N) by recombining the N single skyrmions. Note 
that the sign of this term in the antiferromagnet or- 
dering is negative, therefore smaller seperation between 
skyrmions, R, is favorable. Such a combination costs the 
Coulomb energy. Since the Coulomb interaction is long 
range and the XY term contributes to the total energy as 
a short range interaction, and the order of the Coulomb 
energy and the short range XY interaction with respect 
to the number of particles are 0(N 2 ) and 0{N) respec- 
tively, then the stability of the lattice at the limit of large 
N as a local minima is guaranteed (the collapse of the 
lattice costs too much energy). For large N, an antifer- 
romagnet ordering between the single skyrmions within 
a square lattice can be the minima of the interaction en- 
ergy. Conversely, for the short range interactions, namely 
the hard-core model 



E c [m] oc - J d 2 r(d x m(r) x d y m(r) f , 



(30) 
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we obtain E(Q — N) < NE(Q = 1), i.e., a charge- 
Q skyrmion is .-the global minima, amongst all other 
configurations .E-0 These models are convenient for the 
nuclear physics, to describe the combination of the indi- 
vidual nucleons to make a stable nuclei with the same 
baryonic number. Here the baryonic number is ana- 
logues of the skyrmions charge. This configuration which 
is favorable for the short range strong interaction be- 
tween the nucleons, can be destroyed by the long range 
Coulomb interaction. Therefore the long range nature 
of the Coulomb interaction is crucial for the stability of 
the lattice. Varying the functional form of the inter- 
action term, in E Cl changes the optimized value of the 
self energies for different models of skyrmions and subse- 
quently changes the effective potential among them. A 
cross over between a collapsed phase (no long range or- 
der phase) into a Skyrme lattice occurs due to changing 
the functional form of interaction, from short range to a 
long range. This cross over can be identified by the sign 
of the energy difference (the skyrmions self energy dif- 
ference) between the different topological spin configura- 
tions, E(Q)-E(Q 2 )-E(Qi) where Q = Qi+Q 2 - We end 
up this subsection with a comment. The results of our 
calculation show V e g is negligible as g tends to zero. For 
the zero Zeeman energy the Belavin-Polyakov solutions 
are the well known analytical scale invariant solutions, 
i.e., any analytical function possess the same energy and 
neith.ej-. the s kyrmion positions nor the skyrmion size are 
given. Ej'tSEZl In these circumstances, all of the possible 
Q-sector configurations of a multi-skyrmion system are 
degenerate, e.g., the energy of the N(= Y^iQi) sm gle 
skyrmions is equivalent to the charged- N skyrmion. The 
seperated skyrmions do not feel each other, hence there 
will be no effective potential among them. V e s is non- 
zero if g 7^ 0. 



C. Triangular lattice of the bi-skyrmions 

An estimate based on the total energy of a pair of 
skyrmions with the long range Coulomb interaction, 
shows that the energy of a skyrmion with charge two is 
lower than the energy of two skyrmions, each with charge 



one when g < g c (= 5.3 x 10~ e /e£q).c3 For small values 
of the Zeeman energy a pair of skyrmions is perefered. 
Moreover, a bi-skyrmion (charge-2 skyrmion) is stable to 
accretion of further charge. One may expect a Skyrme 
lattice is governed via the formation of the bi-skyrmions 
when g < g c . The inter-skyrmionic interaction between 
the bi-skyrmions can be obtained by the same calcula- 
tion which has been introduced in Sec. 



methodla to find out the asymptotic behavior of a circular 
symmetric bi-skyrmion. Far from its core we find 



III A. In this 



case, the relation between the hedge-hog fields of the 
bi-skyrmions relative to the standard hedge-hog fields 
can be defined by £j — 2tp — \j- Similar to the case 
of the single-skyrmion lattice, we should have the correct 
shape of a charge-2 skyrmion to obtain the effective po- 
tential. We have solved Eq. (H) , using the finite-difference 



( r ) = Trr-M^r) 



2tt 



(31) 



where Ki{x) = ^Jn/(2x)e~ x (\ + 0(1/ x)) is the modified 
bessel function. Our numerical calculation gives c = 79. 
It is straightforward to show that the interaction between 
the bi-skyrmions can be obtained by: 

2„ [ j2„>Pi( r )Pi( T ') 



VcS { m } = —J2 [d 2 r fd 



r — r 



eg 

' 8tt 3 



J2 cos fe ~ Xi)K (K\Rj - Ri|). (32) 

<ij> 

We observe that the effective Coulomb interaction in 
the multi bi-skyrmionic system is a factor of 4 of the 
single-skyrmionic lattice. Moreover, the sign of the XY 
term is negative, rendering the possibly ferromagnet state 
(Xi = Xj) within the triangular lattice for small values 
of the Zeeman splitting factor. Th is res ult needs more 
comments. We have seen in Sec. [II A that the effec- 



tive potential for the interaction between two seperated 
single-skyrmions is similar to a pair of classical dipolcs, 
i.e., a single-skyrmion behaves as a source of the dipole 
fields. On the other hand, the interaction form of a sys- 
tem of coupled bi-skyrmions acts as a pair of quadrapole 
fields (a charge-2 skyrmion is a source of quadrapole 
fields). This is shown in Fig. ^. The minimum energy 
configuration is therefore ferromagnet alignment of the 
quadrapoles. This can explaine why the negative sign in 
the front of the bi-skyrmionic XY effective potential docs 
exist. However, the rest functional form of the interac- 
tion between the charge-Q skyrmions is Q-independent. 
This stems from the behavior of the asymptotic form of 
our linear field theory,EZI which has been implemented 
to derive the effective interaction among the skyrmions. 
Note that the leading XY-tcrm in |3^ is the quadrapole- 
quadrapole interaction, where the XY dipole-dipole term 
is absent. This is originated from the non-trivial topolog- 
ical structure of the charged-2 skyrmions. The prediction 
of our classical field theory, therefore, rendering a struc- 
tural phase transition, single-skyrmionic square lattice 
into bi-skyrmionic triangular lattice at g c . We expect 
a first order phase transition occurs at g c where a latent 
heat is needed. However, one may expect this phase tran- 
sition can be affected by other mechanisms, like skyrmion 
scattering from the impurities and the quantum disor- 
dering. Since the energy differences needed for observ- 
ing this structural phase transition are very small, due 
to extremely small Zeeman energy, then one expects the 
quantum fluctuations which are present at zero temper- 
ature destroys the long range-order of the bi-skyrmion 
lattice. This requires further investigations, one should 
include such processes to the present formalism in or- 
der to obtain a realistic phase diagram. The results of 



G 



our calculation for the stability of the bi-skyrmion lat- 
tice against the quantum fluctuations is still on the way 
and will be presented elsewhere .E3 Note that the sign of 
the XY-term in Eq.([j2|) in the ferromagnet ordering is 
negative. The collapsed N-charged skyrmion is therefore 
prefered for the short range interaction models, like the 
hard-core model. 



D. Triangular lattice of the single-skyrmions 

Obviously, the skyrmion density may be controlled via 
the Landau level filling factor. For larger separations, 
the XY term in Eq.(p8[) is therefore negligible compared 
to the effective Coulomb energy. This implies that for a 
lower density of the single skyrmions the Coulomb inter- 
action determines the lattice structure and leading to a 
triangular lattice, assuming g > g c . However, for a more 
dense system one may expect a different stable structural 
symmetry, as we have mentioned before. Unlike the case 
where the competition between the monopole Coulomb 
term and the gradient dipole terms is crucial to the sta- 
bility of the square lattice, the electrostatic monopole in- 
teraction is the dominant term to determine the structure 
of the lattice for an extremely dilute [y ~ 1) skyrmions, 
where the Coulomb interaction is long range and the XY 
term contributes to the total energy as a short range in- 
teraction. In this situation a specific skyrmion moves 
in a background which is being made by the vacuum of 
the other skyrmions. Since u t which has been defined 
by Eq.(17), is zero everywhere but in the zth region, 
then p(r) — X^Pi( r )- Clearly the total charge is the 
summation upon the individual skyrmion's charge, i.e., 
Qtot = X)i Qi where Qi is the topological charge of local- 
ized skyrmion in that region. In this case, the interaction 
between separated skyrmions is independent of the rel- 
ative orientation. The skyrmions may be considered by 
point particles where pi(r) = v8(r — R$) and 



Vol 



2c 



d 2 r , pj^wy) 



2 2 N 



r — r 



2c 



^3 



EL — R, 



(33) 



Therefore the ground state is being described by Eq. ( |33] ) 
is clearly a triangular lattice, independent of the 
skyrmion orientation, i.e., a Wigner crystal.EU This is i 
agreement with the field theoretic result of Green et al 
and the microscopic pictures .E2IE3 One may find the crit- 
ical filling factor v\ c where a transition into square lat- 
tice takes place. We define a domain of separation for 
which the square lattice is valid, e.g., vi c < v < v\ c or 
< R < Rq. Here the length scale cut off is denoted 
by Rq, beyond that a phase transition to a triangular lat- 
tice occurs. Note that we have made an approximation 
associated with R > k^ 1 to demonstrate the existence 



of the square lattice. We may expect the approxima- 
tion, namely the linear superposition rule, Eq.(|2C|), fails 
for a high density skyrmions. The crossover between 
low and high density of skyrmions occurs at kR ~ 1 
and leading to the lower critical filling factor i^o i-e., 
V2 C = 1 — g/2ps- A comparison between the ground state 
energy per skyrmion for the square lattice and for the 
triangular lattice yields an estimate for the critical filling 
factor, vi Cl the square-into-triangular transition point at 
<?(> g c ). The result obtained from this calculation leads 
to a condition for the stability of the square lattice 



-3/4 I h 

So I 1 



v\ 1 ^ 4 exp 



-6.333* 



,9o 



> 2 x 10" 



(34) 



Here go = g / (e 2 / elo) ■ Using the typical Zeeman energy, 
g = 0.015e 2 /d?o for GaAs and taking advantage of the 
skyrmion-antiskyrmion duality, we obtain \ v>i c — 1| ~ 0.02 
and \v2c~ 1| ~ 0.3. Our estimate for the structural phase 
transition at v\ c is a factoif-Qf 10 smaller than the micro- 
scopic Hartree-Fock resulted, although it shows qualita- 
tive agreement with that microscopic picture. 



E. Quantum phases 

The appropriate quantum mechanical spectrum of the 
Skyrme lattice in our model can be obtained through 
Eq.(p8|) and the Wess-Zumino action 



N 



X>(|R,-R,I) 

i¥=3 



J2 J(|Ri-Ri|)cos(xi-») U([R;x];r) 



(35) 



where [R; \] = (Ri; xi, R2; X2, — , Rjv; Xn) and R; = 
(Ri cost/?i, Ri siinpi). (pi is the standard azimuthal angle, 
indicating the ith skyrmion. A = ^ 2 £o/(35) is the mo- 
ment of inertia of a single skyrmion, and £ = E z /(2g). 
Obviously, Eq. (|35|) describes a quantum rotor problem 
for x-degree of freedom.lla This equation which has been 
obtained from a first-principle calculation is consistent 
with the result of Ref. 27, and 28 which has been es- 
timated phenomenologically via the fitting the known 
Hartree-Fock long wave length spin waves spectrum. The 
ground state described by Eq. ( |35| ) for skyrmions orienta- 
tion is a quantum antiferromagnet where the U(l) sym- 
metry is broken spontaneously. This model has been used 
to describe the superconductor-insulator (SCI) quantum 
phase transition in panular superconductors and Joseph- 
son junction arrays,t3 and predicts presumably a continu- 
ous quantum phase transition in QHFE3 The appropriate 
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parameters which control the SCI phase transition at zero 
temperature is the Zeeman splitting (g) and the Landau 
level filling factor. Using the empirical data for practical 
systems (g = 0.015e 2 /e£ ~ 2K) shows the XY coupling 
constant, J = (c 2 g/ATr 2 )K (y/ gj2-Kl\p s R), is a fraction 
of the U = {h 2 /2A ) ~ OAK within v > v lc . It leads to 
the prediction of the Mott-insulating phase for this sys- 
tem where we assume that the position of skyrmions are 
fixed, e.g., they are crystalized. Here E z — 0.214e 2 /e^o 
is the Zeeman energy associated with the optimal single 
skyrmion solution corresponding to g = 0.015e 2 /e^oi is 
evaluated by solving Eq.(||), for a single skyrmion. The. 
details of this calculation have been presented elsewhere.EI 
We also expect that the system go through supercon- 
ducting phase for J greater than U since the skyrmions 
earring the electrical charges. In Fig. ||, J and U are 
depicted in terms of the Landau level filling factor for 
a given Zeeman energy. One may see that at v = 1 
the XY-interaction between the skyrmions is negligible. 
As it is shown in Fig. ^ there are two symmetric crit- 
ical point around v = 1. This symmetry is due to the 
duality of the quantum Hall skyrmion-antiskyrmion at 
v = 1. The nature of this phase transition could be 
anticiapted .*ia the bosonic Hubbard model of the quan- 
tum rotors .L3 Its critical behavior can be characterized 
via a map onto a 2 + ID classical XY modcl.t3 Note 
that £o is a continuous variable. It changes continuously 
with the size of skyrmions and/or the Zeeman energy. 
For the half integer value of £o where 2£o is odd, our 
model leads to the Mott-insulating phase about v = 1 
for the quantum Hall liquid where the classical ordered 
phase is destroyed by the quantum fluctuations. Within 
the quantum rotors model, the skyrmions are analogous 
to the vortices. This is a dual picture for the quantum 
Hall state, where the Mott-insulating gap of the bosons 
and the super-current of the vortices are equivalent to 
the quantuia Hall gap and the quantum Hall current, 
repectivelyBa, consistent with the microscopic picture of 
Ref. 27. The energy gap of the Mott-insulating phase, 
A = 2U(= 6g 2 /E z ), is the analogous of the long wave 
length (k — > 0) energy gap of the magnons. Magnons 
are the appropriate collective modes of the quantum Hall 
liquid corresponding to a single flipped spin, propagating 
throught the system. The numerical study for the single 
skyrmionsa shows that E z cx (e 2 /e£o) 2 ^ 3 5 1 ^ 3 as 9 — ► 0. 
This leads to the prediction of vanishing of the Mott- 
insulating gap as A — > g 5 / 3 . We estimate A ~ IK at 
g = 0.015e 2 /e£o- Therefore for the zero Zeeman splitting 
(U = J = 0), the effective Hamiltonian is reduced to the 
exchange Coulomb interaction among the electrons. In 
this case, the gapless Goldstone modes corresponding to 
the spontaneous global SU(2) symmetry breaking (spin 
waves) are responsible for the incompressibility of the 
quantum Hall liquid. The stability of the quantum Hall 
state can be investigated via evaluating the dispersion 
curve of the magnons. 



F. Melting of the Skyrme lattice 

In the following subsection, we implement the theory 
of elasticity^ to study the non-zero temperature phase 
transitions, and the zero-temperature quantum fluctu- 
ations. The universality class of the phases associated 
with the Hamiltonian, Eq.(35), can be changed when 
T > 0. The Skyrme lattice can be anticipated via the 
long-range order of the position of the skyrmions as well 
as the orientational long-range order of their hedge-hog 
fields. These long-range orders can, be destroyed at fi- 
nite temperature. One may show ,£3 the fluctuations of 
the skyrmion orientation fields diverge logarithmically, 
(x 2 ) ~ T\n{L/R), where L is the system size. This char- 
acterize the quasi-long range order of the skyrmions ori- 
entation fields on the 2D-plane. Forming a pair of the 
vortex-antivortex from the hedge-hog fields, leads to a 
finite temperature Kosterlitz-Thouless phase transition 
(Tkt), assuming that skyrmions are crystalized.c3 We 
may predict a critical temperature for the KT phase tran- 
sition based on the classical 2D XY-model. Our calcula- 
tion for the classical 2D nearest neighbor XY-model on 
the square lattice, shows that the phase transition oc- 
curs around T c = 0.9 J ~ AK at v = 0.8 (B ~ IT). 
One may also expect the 2D-Skyrme lattice melts (at 
T = T m ) via the Kosterlitz-lChouless-Halperin-Nelson- 
Young (KTHNY) mechanism.EJ In this case the disorder 
phase can be emerged via mediating a pair of disloca- 
tions. Similarly, the fluctuations of displacement fields 
diverge logarithmicallyEl (R 2 ) ~ T\n(L/R). This ex- 
hibits a finite temperature quasi-long range positional 
order. Our classical field theoretic results for the Tkt 
and T m are shown in Fig. [|. The details of these cal- 
culations will be presented elsewhere.E 2 ] The numerical 
calculation shows that the KT phase transition can not 
account for this system except for small |1 — v\ (v > v\ 
and/or large g where J < /c_BT m , as depicted in Fig 
In contrast to the case of small Zeeman energy where the 
KTHNY is responsible for destroying the translational 
long range order (and therefore melting the Skyrme lat- 
tice) , the orientational long-range order of the hedge-hog 
fields can be destroyed through the KT transition when 
the Zeeman energy is large. However, more accurate cal- 
culation is needed to obtain the melting temperature of 
the Skyrme lattice. For this, we perform a Monte Carlo 
simulation. The results will be presented elsewhere. By 
taking the zero point quantum fluctuations into account, 
the Hamiltonian Eq.(|35|) leads to (R 2 ) ~ |1 — v\a 2 /5 at 
T = 0, where a is the lattice spacing. Using the Linde- 
mann criterion, we observe that the zero point quantum 
fluctuations destroys the long range order of the Skyrme 
lattice when |1 — v\ > 0.5. This is the high density limit 
of the skyrmions where our linear field theory may fail. 
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IV. CONCLUSION 

In this paper we presented the analytical form of the 
inter-skyrmionic interaction via a first-principle calcula- 
tion within the classical field theory. The effective Hamil- 
tonian demonstrated that for extremely dilute skyrmions, 
the ground state is a triangular lattice to minimize the 
Coulomb repulsion. However at higher densities, a square 
lattice forms to optimize the spin gradient and Zee- 
man energies. This work has provided a possibility 
that a triangular to a square lattice phase transition 
which has been previously derived, using the Hartree- 
Fock method, can also be revealed by the field theory 
picture. In contrast to Ref. 11, we show that the field 
theoretic description of the NLcr-model does indeed lead 
to quantitative conclusion in agreement with microscopic 
Hartree-Fock calculations. We also argued the possi- 
bility of the superconducting-insulator phase transition. 
We have shown that within some special short range 
interactions the level-crossing between different kind of 
skyrmions may occur to destablize the Skyrme lattice, 
and a single-skyrmion with the same topological charge 
can be formed. We have shown that for extremely small 
Zeeman energy, skyrmions bound in pair to make a trian- 
gular lattice with the ferromagnets orientation. The pos- 
sibilty of the Kosterlitz-Thoulcss phase transition as well 
as the Kosterlitz-Thouless-Halperin- Nelson- Young phase 
transition has been investigated. 
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FIG. 2. Minima configuration of the system of a coupled 
single- (a) and bi-skyrmions (b) are shown. The single (bi) 
skyrmion act as a source of the classical dipole (quadrapole) 
fields, far from its core. 



FIG. 3. The hopping term of the Bosonic Hubbard model, 
J, and the charging energy U as a function of the Landau level 
filling factor are depicted for a typical Zeeman splitting factor. 



FIG. 4. The prediction of the classical field theory for the 
disordering temperature as a function of the Landau level fill- 
ing factor is shown for g = 0.015e 2 /e£o and g — 0.1e 2 /e£o (the 
inset). In contrast to the case of small Zeeman energy where 
the mechanism of the KTHNY is responsible for the melting 
of the Skyrme lattice, the orientaional long-range order of the 
hedge-hog fields can be destroyed through the KT transition 
for the large Zeeman energies. 
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